Diffusion in a strongly coupled magnetized plasma 
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A first-principle study of diffusion in a strongly coupled one-component plasma in a magnetic field 
B is presented. As in a weakly coupled plasma, the diffusion coefficient perpendicular to the field 
exhibits a Bohm-like l/B-behavior in the strong-field limit, but its overall scaling is substantially 
different. The diffusion coefficient parallel to the field is strongly affected by the field as well and 
also approaches a 1/B scaling, in striking contrast to earlier predictions. 
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The question of how a magnetic field influences the 
transport properties of ensembles of charged particles is 
of considerable significance for a wide spectrum of physi- 
cal systems and situations. Of particular current interest 
are the diffusion properties in the case of strong coupling, 
i.e., when the Coulomb interaction exceeds the kinetic en- 
ergy of the particles. In astrophysics, for example, the 
diffusion of strongly coupled ions directly influences the 
age estimates of white dwarf stars through the timescale 
of gravitational energy release pQ . Also, knowledge of the 
diffusion coefficient of strongly correlated nuclei is a key 
to understand the properties of the crust of magnetized 
neutron stars [2]. Furthermore, the diffusion coefficient 
is directly connected to the energy loss (stopping power) 
of slow ions in dense plasmas [31 H] and, therefore, of spe- 
cial interest for magnetized target fusion scenarios [5] as 
well as in ion beam cooling setups [5J. Finally, recent 
experimental advances with ions in traps as well as with 
dusty plasmas (THSj are starting to access the microscopic 
motion of correlated charge carriers in a strong magnetic 
field. A second line of research originates in the compar- 
atively younger field of strongly coupled plasmas (SCP). 
Here, correlations between particle movements give rise 
to many new effects including non-Fickian diffusion in 
two-dimensional geometries (TUj QT] or shear waves [H] . 
While the mass transport in an unmagnetized SCP has 
been investigated in detail, e.g. Refs. [131 [H], only lit- 
tle is known about transport in plasmas which are both 
strongly coupled and magnetized which is the subject of 
the present paper. 

Recall that diffusion can be understood as a stochastic 
process where particles advance in space on average a 
distance Ar until they undergo a scattering event with 
mean frequency v, giving rise to the diffusion coefficient 

D = (Ar) 2 i/. (1) 

A strong magnetic field B substantially alters diffusion of 
a weakly correlated one-component plasma (OCP). The 
cross-field diffusion coefficient is obtained by replacing 
Ar by the Larmor radius r c ~ B^ 1 and using for v the 
collision frequency. This gives rise to the well-known scal- 
ing of D± ~ B~ 2 [15] which also applies to non-neutral 
confined plasmas if multiple collisions are included [16] . 
A different scaling, D± ~ £? _1 , was predicted first by 



Bohm ("Bohm diffusion") and plays a key role in mag- 
netic fusion plasmas, e.g. Ref. [T7]. More generally, this 
"anomalous" diffusion is expected to be dominant in plas- 
mas with strong field fluctuations, instabilities or plasma 
turbulence, as was shown based on mode coupling theory 
by Marchetti et al. [THj- An overview of the predicted 
scalings is given in Tab. [I] In contrast to weakly cou- 
pled plasmas where diffusion along B is unaffected by 
the magnetic field, in the case of strong correlations par- 
ticle motion perpendicular and parallel to B become cou- 
pled. Suttorp and Cohen have carried out an extensive 
analysis based on the Balescu-Guernsey-Lenard (BGL) 
kinetic equation and predicted that, for strong fields, Dm 
would saturate at two-thirds of the field- free value Do, 
for arbitrary coupling |19) . Finally, there exist molecular 
dynamics (MD) simulation results for a SCP by Bcrnu 
which, however, are limited to three values of B and do 
not allow for a conclusive deduction of scaling laws |20j . 
A more recent MD study [21] reported an unexpected 
non-uniform decay of the diffusion coefficients with B and 
a rapid breakdown of diffusion for large B [2"2] . 

In this work, we report on extensive novel first- 
principle simulations of diffusion of magnetized strongly 
coupled three-dimensional OCPs in a broad range of pa- 
rameters with the following main results (see also Tab. [I): 
1) At large B, a Bohm-like scaling of D± exists at all cou- 
plings, 2) Dm exhibits a Bohm-like scaling as well but only 
in the strongly correlated fluid regime, 3) for moderate 
coupling the scaling of Dm is slower than l/B, however, 
in all cases the decay continues algebraically at strong 
fields and does not saturate as predicted in Ref. |19j . 
These results are of direct relevance for all SCP that are 
well described by the OCP model. 

We now turn to the description of our results. The 
thermodynamic state of the OCP is characterized by the 
coupling parameter T — q 2 x (ATrsgakrjT^ 1 , where a = 

1 /3 

[3/ (4ri7r)] 1 is the Wi gner-Seitz-radius with the number 
density 77, and q and m denote particle charge and mass. 
The strength of the magnetic field B (B || e z ) is given 
by fj — Wc/w,, oc B, the ratio of the cyclotron frequency, 
uj c = qB/m, and the plasma frequency, to 2 — nq 2 / (eqitl). 
In the following, we use reduced units f = r/a and i = 
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TABLE I: Reported B-field scalings of the cross-field and parallel diffusion coefficient of an OCP. Results of this work for 
strong coupling, T > 1, are in bold, see also Tab. |ll| /3 = uj c /uj p , and «fc,6fc,afc,7fc are B-independent coefficients. 



tujp for lengths and time. The equations of motion, 
n = Fi/m + u> c n X e z% i = 1...N, 



(2) 



where Fi is the Coulomb force due to all particles j ' ^ i, 
are solved for N = 8196 particles by standard MD tech- 
niques with Ewald summation and integrators adopted 
to the influence of the magnetic field [231 EHi ■ Periodic 
boundary conditions are imposed in all three directions of 
the cubic simulation box. Prior to measurement, the sys- 
tem is brought into equilibrium by a repeated rescaling 
of the particles' momenta towards the target value of F. 
After that the system is advanced according to Eq. ([2]), 
resulting in a microcanonical propagation. From the 
equilibrium dynamics the diffusion coefficients are cal- 
culated by the Einstein relation for the mean-squared 
displacements [25] 



D\ = lim 
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where (...) denotes averaging over the entire particle en- 
semble 26 . We verified that the results obtained from 
Eqs. ([3j4 ) agree within the statistical uncertainty, with 
the time integral over the velocity auto-correlation func- 
tion Z(t) = (v(t) ■ v(t ))/N d , (Green-Kubo relation [2"5]). 

Consider first the results for the cross-field diffusion 
coefficient D± . Increase of T leads to a rapid monotonic 
reduction of D±, cf. Fig. [T]i. Similarly, a monotonic 
reduction of D± with B is observed, as expected from 
weak-coupling arguments. The agreement with the re- 
sults of Daligault [14] (/3 = 0) and Bernu [20] is excellent 
for B = and fair for magnetized systems. We also in- 
clude kinetic theory results from Ref. [T5] for different 
approximations for the memory kernel, nameley Landau, 
Rostoker and modified Rostoker. While deviations for 
the former two are considerable, the latter turns out to 
be fairly accurate, for B = 0. 

In Fig. [2^, we detail the depencde D±(f3) at different T. 
For L = 1, we observe ^-independence, for small /3, with 
a crossover to a /3 _2 -decay around j3 — 0.1. For large 
/? the decay of D± approaches a asymptotics, as 
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TABLE II: The modulus |a| of the decay exponent of D± and 
D11 (~ P~ a ) averaged for different ranges of j3. The values in 
parantheses correspond to the region [/3 = 5.0 — 100.0]. 



observed for weak coupling in Ref. [TH] , see also [2] . Con- 
sider now the results for strong coupling, T = 2 . . . 100. 
Irrespective of the coupling, D± is independent of f3 for 
small values of (3. With increasing T the intermediate 
f3~ 2 decay is gradually lost. Finally, at high couplings 
(r > 10), a Bohmian regime ~ immediately fol- 
lows for (3 > 0.5. Our findings for the decay exponent a 
(D± r~j P~ a ) are collected in Table[n]for different ranges 
of the magnetic field. Finally, we depict in Fig. [3] the 
ratio of D± to the field-free coefficient Do for ft = 0.5 
and 1.0. At small T, cross-field diffusion is efficiently 
suppressed by the magnetic field. As, however, T in- 
creases, the increased collisionality of the plasma allows 
for a more frequent crossing of the magnetic field lines 
and the diffusion coefficient attains a sizeable fraction of 
its field-free value. 

The present results are, to our knowledge, the first 
observation of Bohmian diffusion in a SCP. Such a be- 
havior is remarkable because it indicates striking simi- 
larities of these highly ordered systems with the com- 
pletely different turbulent fusion plasmas, despite the 
entirely different range of coupling parameters: in both 
cases the dominant transport mechanism are collective 
modes. In an SCP these are predominantly magneto- 
plasmons with a frequency around Co p ~ (uj 2 + cj 2 ) 1 / 2 . 
Indeed, replacing, in Eq. <flj> , Ar — > r c and v uj p yields 
D ± -> j ± (r)k B T(qB)- 1 lor (3 > 1. For example, for 
T = 100 we find 7j_ ~ 2/3 which is surprisingly close to 
Spitzer's weak coupling estimate of 7 = 1/2.2 [T7] but 
substantially larger than the familiar 1/16 of Bohm. 
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FIG. 1: (color online) Diffusion coefficient a) perpendicular and b) parallel to the magnetic field as a function of coupling for the 
unmagnetized system and two moderate values of B. Also shown are the results of Daligault [T^j (solid black line), Bernu [20) 
(crosses) and from kinetic theory (BGL [19]). 




FIG. 2: (color online) Diffusion coefficient a) perpendicular and b) parallel to the magnetic field versus field strength ft = w c /aj p 
for different F. The straight broken lines indicate a decay of f3~ 2 , /3 _1 , and /3 -0 ' 3 , respectively. Also shown are simulation 
data from Ref. [2T] for T = 50 (open diamonds) [22] • The open circles in b) are predictions from BGL theory with a modified 
Rostoker Kernel for F = 1 and 10 [19]. 



Let us now turn to the field-parallel diffusion. As ex- 
pected, for small couplings, Dn is practically insensitive 
to B, (Fig. [TJd). In contrast, in an SCP in the liquid 
state that sustains shear, D\\ decreases with increasing 
P. This effect is elucidated in more detail in Fig. ^show- 
ing the ratio Dy/Do as a function of T for two values 
of j3. At small T, the ratio is close to unity and de- 
clines steadily with growing T until, around r crit w 30, it 
becomes T-independent [in Fig. [2jb this corresponds to 



parallel curves for different T]. Interestingly, the onset 
of this saturation coincides with a change in the particle 
dynamics: at r cr it, the velocity autocorrelation function 
Z(t) [see inset of Fig. [4] changes from a modulated decay 
to large amplitude oscillations which extend to negative 
values [we demonstrate this by showing, in Fig. |4j the 
integral over the negative portions of Z(t) from t = to 
t = 100]. Such negative correlations reflect "caging" of 
particles in their local potential minima which lasts on 
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FIG. 3: (color online) D± in units of its field-free value (for 
the same F) for two field strengths, /3 = 0.5 and 1.0. 
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FIG. 4: (color online) D\\ in units of its field-free value (for 
the same F) for f3 — 0.5 and 1.0. Also shown are the in- 
tegrated negative velocity auto-correlations (caging time, see 
text). Inset: The velocity autocorrelation function Z(t) for 
three values of T. 



average a time T c . Donko et al. showed that, in this 
regime, D ~ T^ 1 [27 . Thus, Fig. indicates that m an 
external B field parallel diffusion for Y > Y clit is limited 
by caging effects whre the magnetic field prolongs T c . 

The /3-dependence of D\\ is detailed in Fig. j2jj and 
Tab. |n] In all cases, an initial /3°-behavior is followed by 
an algebraic decay, Dm ~ /3~ Q . The exponent a is smaller 
than unity for comparatively small Y. With increasing T, 
however, a gradually increases and reaches a ~ 1 at Y sa 
r cr it(/3)- Thus, there exists Bohmian diffusion parallel to 
the field which is another novel feature of SCPs. Note 
that this rules out the predicted saturation D\\(/3)/Do — » 
2/3, for large (3 of Ref. [19]. 

In conclusion, we have presented first-principle simu- 
lation results which, for the first time, illuminate the dif- 
fusion characteristics of SCPs under external magnetic 
fields in great detail, substantially generalizing results 
known for weakly coupled OCPs. At low coupling, three 
regimes of the magnetic-field dependence of D± have 
been identified, with a faster decay at intermediate B- 
fields. At large coupling parameters and strong magnetic 
fields we have found evidence of Bohmian 1 / fl-diffusion 
indicating the dominant role of collective modes. On the 
other hand, field-parallel diffusion in a SCP differs drasti- 
cally from a weakly coupled plasma: Dy also depends on 
B, resulting in an algebraic decay, Dm ~ B~ a . This is a 
unique feature of high-T liquid-like OCPs that is induced 
by many-particle correlations - in particular shear and 
caging effects. Our results are expected to be of direct 
relevance for the SCP in compact stars and in inertial 
confinement fusion as well as for trapped ions and dusty 
plasmas. Moreover, they suggest to look for similar be- 
havior of other transport quantities of magnetized SCP. 
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